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Abstract: The antenna subtraction method for the computation of higher order correc- 
tions to jet observables and exclusive cross sections at collider experiments is extended to 
include hadronic initial states. In addition to the already known antenna subtraction with 
both radiators in the final state (final- final antennae), we introduce antenna subtractions 
with one or two radiators in the initial state (initial- final or initial-initial antennae). 
For those, we derive the phase space factorization and discuss the allowed phase space 
mappings at NLO and NNLO. We present integrated forms for all antenna functions 
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terms at NLO on two examples. The extension of the formalism to NNLO is outlined. 
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1. Introduction 



The calculation of perturbative higher order corrections to exclusive observables (especially 
jet production cross sections, but also transverse momentum or rapidity distributions) re- 
quires a systematic procedure to extract infrared singularities from real radiation contribu- 
tions. These singularities arise if one or more final state particles become soft or collinear. 
For the task of next-to-leading order (NLO) calculations [1] several systematic and process- 
independent procedures are available. Except for the phase space slicing technique [2], all 
these methods [3-7] work by introducing additional terms, that are subtracted from the 
real radiation matrix element at each phase space point. These subtraction terms approx- 
imate the matrix element in all singular limits, and are sufficiently simple to be integrated 
over part of the phase space analytically. After this integration, infrared divergences of the 
subtraction terms are made explicit, and the integrated subtraction terms can be added 
to the virtual corrections, thus yielding an infrared finite result. One of these subtraction 
methods is antenna subtraction [5,6], which constructs the subtraction terms from so-called 
antenna functions. These antenna functions describe all unresolved partonic radiation (soft 
and collinear) between a hard pair of radiator partons. 

Extensions to next-to-next-to-leading order (NNLO) are discussed in the literature for 
phase space slicing [8], and for several subtraction methods [9-14]. A completely inde- 
pendent approach, avoiding the need for analytical integration is the sector decomposition 
method, which has been derived for virtual [15] and real radiation [16] corrections to 
NNLO, and applied to several observables already [17]. Among the subtraction methods, 
only the NNLO formulation of the antenna subtraction method [14] is worked out to a 
sufficient extent to be readily implemented in the calculation of NNLO corrections to a 
physical process. Using this method, the calculation of e + e~ — > 3 jets at NNLO accuracy 
is currently under way [18]. However, up to now, the antenna subtraction (both at NLO 
and NNLO) method was developed only to handle unresolved singular radiation off final 
state partons. An extension to radiation off initial state partons was missing so far for this 
method, while all other NLO subtraction methods could handle radiation off initial and 
final state particles. 

It is the purpose of this paper to extend the antenna subtraction method to include 
radiation off initial state partons, such that it can be used in the calculation of higher 
order corrections to processes at lepton-hadron or hadron-hadron colliders. In Section 2, we 
recall the structure of antenna subtraction terms and describe their systematic derivation. 
Antenna subtraction for radiation off final state partons was formulated already in [5,6,14] 
and is briefly summarized in Section 3. In Sections 4 and 5, we present a detailed derivation 
of antenna subtraction for situations with one or two radiator partons in the initial state. 
In these sections, we provide all necessary ingredients for an implementation at NLO, and 
discuss the extension to NNLO. To illustrate the application method, we describe two NLO 
examples in Section 6. Finally, Section 7 contains conclusions and an outlook. 
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2. Antenna subtraction 



To obtain the perturbative corrections to a jet observable at a given order, all partonic 
multiplicity channels contributing to that order have to be summed. In general, each 
partonic channel contains both ultraviolet and infrared (soft and collinear) singularities. 
The ultraviolet poles are removed by renormalisation in each channel. Collinear poles 
originating from radiation off incoming partons are an inherent feature of the incoming 
partons, and are cancelled by redefinition (mass factorization) of the parton distributions. 
The remaining soft and collinear poles cancel among each other only when all partonic 
channels are summed over. 

While infrared singularities from purely virtual corrections are obtained immediately 
after integration over the loop momenta, their extraction is more involved for real emission 
(or mixed real-virtual) contributions. Here, the infrared singularities only become explicit 
after integrating the real radiation matrix elements over the phase space appropriate to the 
jet observable under consideration. In general, this integration involves the (often iterative) 
definition of the jet observable, such that an analytic integration is not feasible (and also 
not appropriate). Instead, one would like to have a flexible method that can be easily 
adapted to different jet observables or jet definitions. Therefore, the infrared singularities 
of the real radiation contributions should be extracted using infrared subtraction terms. 
The crucial points that all subtraction terms must satisfy are that (a) they approximate 
the full real radiation matrix elements in all singular limits and (b) are still sufficiently 
simple to be integrated analytically over a section of phase space that encompasses all 
regions corresponding to singular configurations. 

For NLO calculations, several different methods are available to derive subtraction 
terms in a process-independent way [3-7]. One of these methods is the so-called antenna 
subtraction [5,6]. 

In this method, antenna functions describe the colour-ordered radiation of unresolved 
partons between a pair of hard (radiator) partons. All antenna functions at NLO and 
NNLO can be derived systematically from physical matrix elements, as shown in [14,19]. 
They can be integrated over the factorized antenna phase space [9] using loop integral 
reduction techniques extended to phase space integrals [21], and then combined with virtual 
corrections to partonic processes with lower multiplicity. 

Up to now, antenna subtraction has been formulated at NLO [5, 6] and NNLO [14] 
only for processes with a colourless initial state. In this case, both radiator partons are 
in the final state, we call this situation a final-final antenna. For collider observables 
involving hadronic initial states, there can be either one or both partons in the initial 
state. Unresolved radiation off these initial state partons can also be subtracted using 
antenna functions, with one or two radiators in the initial state. We call these initial-final 
and initial-initial antennae. The radiated parton is always in the final state. Figures 1-3 
illustrate how a single unresolved parton can be emitted between radiators in the final or 
initial state, and show how all these situations are factorized into antenna functions. Each 
antenna contains both collinear limits of the unresolved parton with either radiator as well 
as the soft limit. 
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Figure 1: Antenna factorization for the final-final situation. 
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Figure 2: Antenna factorization for the initial-final situation. 



j 




Figure 3: Antenna factorization for the initial- initial situation. 

In each situation, the subtraction term is constructed from products of antenna func- 
tions with reduced matrix elements (with fewer final state partons than the original matrix 
element), and integrated over a phase space which is factorized into an antenna phase space 
(involving all unresolved partons and the two radiators) multiplied with a reduced phase 
space (where the momenta of radiators and unresolved radiation are replaced by two rede- 
fined momenta). These redefined momenta can be in the initial state, if the corresponding 
radiator momenta were in the initial state. The full subtraction term is obtained by sum- 
ming over all antennae required for the problem under consideration. In the most general 
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case (two partons in the initial state, and two or more hard partons in the final state), this 
sum includes final-final, initial-final and initial-initial antennae. 
To specify the notation, we consider the hadronic cross section 

dcJ = E / ^T^rfa/ii^fb^da^H,,^) , (2.1) 

where £1 and £2 are the momentum fractions of the partons of species a and b in both 
incoming hadrons, with / being the corresponding parton distribution functions and H\ 2 
denoting the incoming hadron momenta. The cases of only one or no incoming hadrons 
are obtained trivially by replacing the relevant /(£) by (5(1 — £). The dependence of the 
parton- level cross section da on the parton species a, b is obvious, and not stated explicitly 
for ease of notation. It should be noted that the parton level cross section da is normalized 
to the hadron-hadron flux factor, which is transformed into the parton-parton flux factor 
by dividing out £1 and £2 in the above. 

We define the partonic tree- level n-parton contribution to the m-jet cross section (for 
tree- level cross sections n = m; we leave n 7^ m for later reference) in d dimensions by, 



n't Pi: 



P2) 



x}|K(fci, • • • , k n ; Pl ,p 2 )\ 2 . . . , k n ). (2.2) 

The definition of the observable Jm\k\, . . . , k n ) can depend on (for example through 
cuts on the jet rapidities), although this is not stated explicitly here. The normalization 
factor N includes all QCD-independent factors as well as the dependence on the renormal- 
ized QCD coupling constant a s , ^ n denotes the sum over all configurations with n partons, 
d<I> n is the phase space for an n-parton final state with total four-momentum p^ + p% in 
d = 4 — 2e space-time dimensions, 

(!*„(*!, • • • , k n ; Pl , P2 ) = 2 ^ i( J;L • • • 2E n {2^ ^ 5d{Pl +P2-ki-.---k n ), 

(2.3) 

while S n is a symmetry factor for identical partons in the final state. |-M„| 2 denotes a 
squared, colour-ordered tree-level n-parton matrix element. 

Antenna subtraction terms da are constructed using parton-level antenna subtraction 
terms as in (2.1), such that 

da - da 5 

is finite in all unresolved limits, and that phase space integrals contained in it can be carried 
out numerically. 

In the following, we will briefly summarize the features of antenna subtraction in the 
final-final case, and derive the antenna phase spaces and antenna functions for the two 
other cases. 
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3. Final-final configurations 

In configurations involving final-final antennae, both radiators are in the final state. This 
case was described previously in detail at NLO [5,6] and NNLO [14]. The NLO subtraction 
term for an unresolved parton j, emitted between hard final-state radiators i and k is 
depicted in Figure 1. It reads 



da s ' {ff \p 1 ,p 2 ) =M ^2d& m+1 (k 1 ,...,k m+1 ;p 1 ,p2)-z 

m+l ^ m+1 

x E X ijk \M m (k!, . . . ,K T , K K , . . . , k m+1 ; Pl ,p 2 )\ 2 J ( ™\ki, ■ ■ ■ , K z , K K , . . . , k m+1 ) 



The subtraction term involves the m-parton amplitude depending only on the redefined on- 
shell momenta k±, . . . , Kj, Kk, ■ ■ ■ , k m+ \ where Ki, Kk are linear combinations of ki, kj, k^ 
while the tree antenna function X^ k depends only on fcj, kj, kk- X^ k describes all of the 
configurations (for this colour-ordered amplitude) where parton j is unresolved. 

The jet function in (3.1) does not depend on the individual momenta ki, kj and 
k/., but only on Kj,Kk- One can therefore carry out the integration over the unresolved 
dipole phase space appropriate to ki, kj and kk analytically, exploiting the factorization of 
the phase space, 

d$ m+ i(fci, . . . , k m+1 ;p 1 ,p 2 ) = 
d$ m (fci, ...,Kj,K K ,..., k m+1 ;p 1 ,p 2 ) • d$x ijk (ki, kj, k k ; Kj + K K , 0) . (3.2) 

The NLO antenna phase space d<&x ijk is proportional to the three-particle phase space 
relevant to a 1 — > 3 decay. 

At NNLO, one has to consider the emission of one parton in a one- loop corrected 
process, or the emission of two partons at tree level. Both these were described in detail 
in [14]. While the one-loop antenna subtraction is largely an extension of the above with the 
replacement of the tree-level antenna function by a one-loop antenna function, X° — ► X 1 , 
several new features appear in the subtraction of two unresolved partons at tree-level. 

In particular, one must pay attention to the colour-connection of the two unresolved 
partons. If they are colour-unconnected or almost colour-unconnected (sharing a common 
radiator), the subtraction term is obtained by iterating the procedure employed at NLO, 
now yielding products of two antenna functions. If both unresolved partons j, k are colour- 
connected, new four-parton antenna functions Xijki appear in the subtraction terms: 



where the sum runs over all colour- adjacent pairs j, k and implies the appropriate selection 
of hard momenta As before, the subtraction term involves the m-parton amplitude 



(3.1) 




x \M m {h, . . . , K l7 K L , . . . , k m+2 ; Pl , P2 )\ 2 (h, . . . , K L , . . . ,k m+2 ) 
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evaluated with on-shell momenta k\, . . . , Kj, Kl, . . . , fc m +2 where now Kj and Kl are a 
linear combination of ki, kj, kk and k\. As for the NLO antenna of the previous section, 
the tree antenna function Xf- kl depends only on ki, kj, kk,h. Particles i and / play the role 
of the radiators while j and k are the radiated partons. 

Once again, the jet function jffl in the above equation depends only on the parent 
momenta Kj, Kl and not hi, . . . ,k[. One can therefore carry out the integration over the 
unresolved antenna phase space (or part thereof) analytically, exploiting the factorization 
of the phase space, 

d$ m+2 (/ci, . . • , k m+2 ;pi,p 2 ) = 
d$ m (/ci, . . . , Ki, K L , . . . , k m +i;pi,P2) ■ d®x ljkl (ki, kj,k k , h ; Ki + Kl , 0) . (3.3) 

This phase space factorization must be carried out such that all unresolved limits are 
reproduced correctly. The most general parameterization for this case is derived in [9]. 
It should be noted that &<&x ijkl is proportional to the 1 — > 4 parton phase space; the 
analytical integration of the antenna functions over this phase space can thus be carried 
out with standard methods [20,21]. 



4. Initial-final configurations 

In the presence of hadrons in the initial state, matrix elements exhibit singularities that 
are not accounted by the subtraction terms discussed in the previous section. These singu- 
larities are due to soft or collinear radiation within an antenna where one or the two hard 
partons are in the initial state. 

As discussed in [14], the terms necessary to subtract singularities associated with col- 
ored particles in the initial state can be simply obtained by crossing the corresponding 
antennae for final state singularities. Due to the different kinematics involved, the factor- 
ization of phase space is slightly more involved and the corresponding phase space mappings 
are different. To cancel explicit infrared poles in virtual contributions and in terms arising 
from parton distribution mass factorization, the crossed antennae must be integrated, an- 
alytically, over the corresponding phase space. In this section we will present the antennae 
and phase space mappings to subtract singularities when only one of the radiating partons 
is in the initial state. 



4.1 Subtraction terms for initial-final singularities 

Subtraction terms in the case of one hard parton in the initial state are built in the same 
fashion as for the final-final case (formula (2.5) in [14]). We have the following subtraction 
term associated to a hard radiator parton i with momentum p in the initial state: 

da S '( if \p, r) = M V d$ m+ i(/ci, . . . , k m+l ;p, r) — — 
m+i bm+1 

x Y. X hk \M m (ki, ...,K K ,..., k m+1 ;xp, r)\ 2 J^Hh, ...,K K ,..., k m+1 ) . (4.1) 
j 
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The additional momentum r stands for the momentum of the second incoming particle, 
for example, a virtual boson in DIS, or a second incoming parton in a hadronic collision 
process. This contribution has to be appropriately convoluted with the parton distribution 
function /j. The tree antenna X^- k , depending only on the original momenta p, kj and 
fcfc, contains all the configurations in which parton j becomes unresolved. The m-parton 
amplitude depends only on redefined on-shell momenta k\ , . . . , Kk , • • • , and on the mo- 
mentum fraction x. In the case where the second incoming particle is a parton, there 
is an additional convolution with the parton distribution of parton r and corresponding 
subtraction terms associated with it. 

The jet function, jffi, in (4.1) depends on the momenta kj and kk only through Kk- 
Thus, provided a suitable factorization of the phase space, one can perform the integration 
of the antennae analytically. Due to the hard particle in the initial state, the factorization 
of phase space is not as straightforward as for final-final antennae. We start from the 
(m + l)-particle phase space 

d<S> m+1 (k 1 ,...,k m+1]P ,r) = (2ir) d 5^r+p-^k)j ]][dki] (4.2) 

where [dfc] = d d fc<$+(ife 2 )/(27r) d - 1 . We insert 

1 = J d d q5(q + p-k 3 -k k ) , (4.3) 

and 

1 = % I T / [AKk] {27T)d6 { - q + xp ~ Kk ^ ' ( 44 ) 

with Q 2 = —q 2 . Finally, integrating over q, the phase space can be factorized in an m- 
parton phase space convoluted with a two particle phase space: 

d$ m+ i(fci, . . . , k m+1 ;p,r) = d$ m (fci, . . . ,K K , . . . ,k m+1 ;xp,r) 

Q 2 dx 
x — d<5> 2 (kj,k k ;P,q)— ■ (4.5) 

Replacing the phase space in (4.1), we can explicitly carry out the integration of the antenna 
factors over the two particle phase space. When combining the integrated subtraction terms 
with virtual contributions and mass factorization terms, it turns out to be convenient to 
normalize the integrated antennae as follows 

■y 



" t = W)I d ^ x -'- (4 - 6) 



where 

p-nE 

C(e) = (4nY . (4.7) 

The integrated form of the subtraction term is then 

da s >W)(p,r) = J^J^^ / —C(e)X i:jk (x)d^ m (k 1 ,...,K K ,...,k m+1 ;xp,r) 

m+l j m+1 J X 

x \M m (h, . . . , K K , . . . , k m+1 ;xp, r)\ 2 J^{k u ...,K K ,..., fc m+ i)(4.8) 
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Finally, the subtraction term has to be convoluted with the parton distribution functions to 
give the corresponding contribution to the hadronic cross section. The explicit poles in the 
integrated form cancel the corresponding ones in the virtual and PDFs mass factorization 
contributions. To carry out the explicit cancellation of poles, it is convenient to recast, 
by a simple change of variables, the integrated subtraction term, once convoluted with the 
PDFs, in the following form 

xC(e) X iJk (x) d<7 B (£iifi,&# 2 ) . (4.9) 

This convolution has already the appropriate structure and mass factorization can be car- 
ried out explicitly leaving a finite contribution. The remaining phase space integration, 
implicit in the Born cross section, da B , and the convolutions can be safely done numeri- 
cally. When considering reactions with only one incoming hadron, the second PDF has to 
be replaced by a Dirac delta. Reactions with two hadrons will require additional subtrac- 
tions containing initial-final antennae involving the second parton in the initial state and 
initial-initial antennae as well. This case will be discussed in Section 5 below. 

4.2 Phase space mapping 

Proper subtraction of infrared singularities requires that the momenta mapping satisfies 

xp — > p , Kk — ► kj, when j becomes soft , 

xp — > p , Kk — ► kj + kk when j becomes collinear with k , 

xp — ► p — kj , Kk — ► kk when j becomes collinear with i . (4.10) 

In this way, infrared singularities are subtracted locally, except for angular correlations, 
before convoluting with the parton distributions. That is, matrix elements and subtraction 
terms are convoluted together with PDFs. In addition, the redefined momentum, Kk, 
must be on shell and momentum must be conserved, p — kj — kk = xp — Kk, for the phase 
space to factorize as above. 

As discussed, in the case of configurations with two hard radiators in the final state, the 
three-to-two-parton map of [9] is suitable, as it treats both collinear limits symmetrically 
and there is only one mapping describing all the singular configurations contained in the 
antennae. 

When subtracting initial state singularities, however, the mapping of [9] leads to a 
non factorizing phase space. The decisive point is that this mapping, modified to account 
for a particle in the initial state, introduces a new initial state momentum, P as a linear 
combination of p, kj and kk- However, as there is no integration over P, the factorization 
of phase space would not be complete, because the m-parton matrix elements depend on 
P. On the other hand, if P is proportional to p, factorization is granted, in the form 
of a convolution between the reduced matrix elements and the integrated antennae, as 
we detailed above. In this case, we inmediately obtain the dipole momentum mappings 
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of [3], combined into a single mapping interpolating between all the singular limits of the 
antennae. Explicitly: 

_ Slj + Sik — Sjk 

x — ) 

Slj + Sik 

K K = kj + k k - (1 - x)p , (4.11) 

where sij = (p — kj) 2 , etc. If parton j becomes soft or collinear to parton k, x — ► 1. If 
parton j becomes collinear with the initial state parton i, x = 1 — z with z the fraction of 
p carried by j. 

The mapping in eq. (4.11) is, in addition, easily generalized to deal with more than 
one parton becoming unresolved. The building blocks for the double real radiation in 
the initial-final situation are colour-ordered four-parton antenna functions Xij^i, with one 
radiator parton i (with momentum p) in the initial state, two unresolved partons j, k and 
one radiator parton I in the final state. Starting with the generalization of (4.3) to three 
particles in the final state, and combining with (4.4) we have the following mapping at 
NNLO: 

Slj + Sik + s ll — Sjk — Sjl — Ski 

x = — , 

■Slj + Sik + su 

K L = kj + k k + h-{\- x)p , (4.12) 

where kj, kk and ki are the three final state momenta involved in the subtraction term. This 
mapping can be obtained from the tripole mapping [21,22] for final-final configurations at 
NNLO. It satisfies the appropriate limits in all double singular configurations: 

1. j and k soft: x — > 1, Kl — > h, 

2. j soft and kk \\ kf. x — > 1, Kl — > kk + h, 

3. kj = zp || p and k soft: x — > 1 — z, Kl — > h, 

4. kj = zp || p and kk || kf. x — > 1 — z, Kl — ► kk + h, 

5. fcj || kk || /fy: x — > 1, Kl — > Aij + fefc + ^, 

6. fcj + fcfe = || p: x — > 1 — z, i^L — > ^, 

where partons j and A; can be interchanged in all the cases. 

The construction of NNLO antenna subtraction terms requires moreover that all single 
unresolved limits of the four-parton antenna function Xijki have to be subtracted, such 
that the resulting subtraction term is active only in its double unresolved limits. A system- 
atic subtraction of these single unresolved limits by products of two three-parton antenna 
functions can be performed only if the NNLO phase space mapping turns into an NLO 
phase space mapping in its single unresolved limits [9]. 

In the limits where parton j becomes unresolved, we denote the parameters of the 
reduced NLO phase space mapping (4.11) by x' and K' L . We find for (4.12): 
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1. j becomes soft: 



sik + su- s k i , v i.i ( i x vi 

■ = x , K L -> fc fc + ki - (1 - x)p = K L . 

sik + su 



2. fcj || fe fc , kj + k k = K K : 
si, 

x — > 

six + Su 



x -> ^ + Su — s JSi = x' , K L — > k K + fe/ - (1 - ar)p = . 



3. fcj = || p: 

X ^ {1 - Z)(sik + Sll) - Skl =(l-z)x', K L ^k k + ki-(l-x')(l-z)p = K' L . 
Sik + su 

It can be seen that in the first two limits, the NLO mapping involves the original incoming 
momentum p, while in the last limit (initial state collinear emission), it involves the rescaled 
incoming momentum (1 — z)p. To subtract all three single unresolved limits of parton j 
between emitter partons % and k from Xijki, on e needs to subtract from it the product 
of two three-parton antenna functions X^jk • Xj ^i- The phase space mapping relevant to 
these terms is the iteration of two NLO phase space mappings. Analytical integration of 
these terms with this mapping will result in a double convolution of both antenna functions 
with the reduced matrix element. 

Equally, parton k can become unresolved. Expressing the reduced NLO phase space 
mapping by x" and K'[. We find for (4.12): 



1. k becomes soft: 



* li + 8 " Sjl =*", K L ^k j + h-{l-x)p = K'[. 
Slj + su 



2. k k || kj, kj + k k = K K : 

sik + su - ski 
x — ^ 

SlK + Su 

3. k k || ki, ki + k k = K K : 



x". 



K L ^k K + k-(l- x)p = K'l . 



x SlK + Slj SK] = x" , K L k K + kj - (1 - x)p = K'l . 

SlK + Slj 

In all limits, the reduced NLO mapping involves the original incoming momentum p. Con- 
sequently, the three single unresolved limits of parton k between emitter partons j and I 
can be subtracted from Xj jkl by a product of a final-final and an initial-final three-parton 
antenna function Xj k \ ■ XijL. The phase space mapping relevant to these terms is the 
product of an NLO final-final phase space mapping with an initial-final mapping. Inte- 
gration of the final-final antenna phase space yields an constant, not involving an extra 
convolution, such that these terms appear in the integrated subtraction term only with a 
single convolution with the reduced matrix element. 
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4.3 NLO antenna functions 

We now present explicit results for all the antenna functions necessary to subtract infrared 
singularities associated with one particle in the initial state. The unintegrated form of 
all of them can be obtained from the corresponding expressions for the tree level three 
particles antennae in [14] by appropriate crossing of particles from the final to the initial 
state. In the cases where there are different particles in the final state, there are more than 
one possible crossing and, thus, more than one corresponding antenna. 

The invariants for antenna X^j^ are defined as s = (kj + A^) 2 , t = (p — kj) 2 , 
u = (p — kk) 2 and Q 2 = —q 2 , where q = p — kj — k^. For the integrated antennae we define 
x = Q 2 /(2p ■ q). The colour-ordered splitting kernels are given by 



V { S{x) = -S(l-x)+2V (x)-l-x, 
p q °J(x) = l-2x + 2x 2 , 



p<$(x) = - -2 + x, 



X 

11 



pf g \x) = y 5(1 - x) + 2V (x) + 3 - 4 + 2x - 2x 2 , 
P ( ggA x ) = -o 8(1 -x) 



(4.13) 



where we introduced the distributions 



V n (x) 



ln n (l - x) 
1-x 



+ 



The colour-ordered infrared singularity operators are as in [14]: 



id) 

<?<? 



( e ) s qq) 



=«7 



1(1) ( 



qg V e ) s qg 



I (1) fe s 
L gg v e ' 6 . 

r(l) 



.'/.'/ i 



qq,F 
j(i) 

^gg,F( e i s gg) 



2T(1 - e) 
2Y(1 -e) [e 



+ 



2T(1 - e) 



3 

2? 
5_ 
37 
11 

6e 



qq 



M(-s qg ) 

K(Sgg) 



(e, S qq ) = 

( e ' s qg) = 



=^7 



2r(i 



i 

6e 
1 

2r(l - e) 3e 



on 



M(- Sqg y 
M(- Sgg y 



(4.14) 



Although the antenna functions are obtained by simple crossing of the antenna functions for 
the final-final case, there are some important differences in the decomposition of antenna 
functions into sub-antennae. In the final-final case, antenna functions involving a hard 
gluon radiating unresolved gluons had to be split into different configurations since any 
final state gluon could be identified as the hard radiator. This ambiguity is no longer 
present if a gluon is crossed into the initial state, since an initial state gluon is hard by 
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kinematical constraints. Instead, a different ambiguity appears, since the initial state gluon 
can split either into a quark or into a gluon, thus leading to two possible reduced matrix 
elements. This ambiguity requires decomposition of the relevant gluon-initiated antenna 
functions into sub-antennae according to criteria completely different from the final-final 
situation, as will be discussed in Section 4.3.2 below. 



4.3.1 Quark initiated antennae 

We consider first antennae with a quark in the initial state. There is one quark-quark 
antenna, given by 



A = -■ 



1 (2u 2u 2u 2 t 



, + T + ir + ; + ii + 0(e) - 



(4.15) 



Its integral over the phase space (normalized as in eq. (4.6)) gives 



A 



2l\ L >(Q 2 )5(l-x) + (Q 2 



3 , . , . 3 x 1 ~\~ x _ , . 1 ~\~ x 

- -V (x) +V 1 (x)-— — log(l -x)- ^—^ log(x) + 0(e) 



4 " v ' y ' 2 
There are three quark-gluon antennae given by 



(4.16) 



D 



q,gg 



s 2 s 2 t 2 At 2 Au 2 Au 2 3st 3su 
— + — + — + — + + — + — + — 

t u s u s tut 



{Q 2 f 

2t 3 2u 3 6tu 

+ — + + + 6s + 9t + 9u + 0(e) , 

su st s 



1 



EX 



(Q 2 ) 2 
l 



t 2 u 2 
— + — + t + u\ +0(e), 
s s ■ 

s 2 u 2 

— + -J + s + u ) +°( e ) 



When integrated over the factorized phase space, they yield 
-4I«(Q 2 )5(l-x) + (Q 2 )- 
1 



(4.17) 
(4.18) 
(4.19) 



9,99 



67 



11 



- -7T 2 8(1 -x) - —V Q (x) 
18 3 / v ; 6 v ' 



+2V i (.(•) - — + 1 - x - (1 + x) log(l - x) - IdJL iog( x ) + 0(e) 
6x 1 — x 



c q,q'q> 



9,99 



4I« (Q 2 )5(l - x) + (Q 2 r 



(Q 



qg 

2\-e 



.|*(i-, ) + i lh (,)-i + ( e ) 



2e 



x 



2x 



1-x 



X 



+ 0(e) 



Finally, there is one gluon-gluon antenna with a quark in the initial state: 

1 



q,qg 



(Q 2 ) 2 



yielding 

Glqg = (Q'r 
when integrated over the antenna phase space 



7 1 (2-2x + x 2 ), 
+ 2» '° 8 



1 — x 



+ 0(e) 



(4.20) 
(4.21) 
(4.22) 

(4.23) 

(4.24) 
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4.3.2 Gluon initiated antennae 

For the gluon initiated antennae, we find one quark-quark antenna 



A 



--^1t + z 7 + i7 + i + L J +0(() - 



2s 2s 2s 2 



t 



Its integrated form is 



A - 
•^am 



(Q 



2\-e 



-pg)(x) - (1 - 2x + 2x 2 )(log(l - x) - log(x)) + 0(e) 



(4.25) 



(4.26) 



There is one quark-gluon antenna with a gluon in the initial state 

1 /„,2 „,2 ^2 ^2 ^„2 ,i„2 
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(Q 2 ) 2 

2t 3 2s 3 



W u" 4t z 4t z 4s z As z 3tu 3su 

— + — + — + — + + — + + — 

t s s u u t s t 



Qst 



+ — + + — + 6n + 9t + 9s)+ 0(e) . 

su tu u 



(4.27) 



This antenna contains singular limits when the quark or the gluon in the final state become 
collinear with the initial state gluon. In the first case it collapses into a quark-gluon antenna 
and in the second case into a gluon-quark one. Accordingly, the reduced matrix elements 
accompanying these two singular configurations are different. Thus, the antenna must be 
split to separate these two configurations. This can be easily done by partial fractioning 
in the variables t and u, we obtain 

1 I /0„,2 Qc ,2 



D 



am 



2(Q 2 ) 2 



2u 2 8s 2 6su 4s 3 

~T + ~T + ~T ~ t{Q 2 + 7) 



+ 0(e), 



(4.28) 



and 



D 



1 1 



am 



2 (Q 2 ) 2 
I2su 



4s 3 



2t 2 8u 2 8u 2 8s 2 6tu Au z 
s s t t s st (Q 2 + s)t 



+ 



+ 12t + 18u + 18s + 0(e) 



(4.29) 



where we have adjusted the names of the antennae so that D®- now does not contain 
singularities when j becomes collinear with the initial state gluon. We also changed the 
sign of the first sub-antenna and exchanged t and u in the second case to agree with the 
definitions given at the beginning of the section. The two sub-antennae can be integrated 
over the factorized phase space, namely: 



V 



am 



(Q 



2\-£ 



-y/S^ + h~ l + \ {l ~ 2x + 2x2) log(1 



X 



- (1 - 2x + 2x 2 ) log(x) + 0(e) 



(4.30) 



and 



Ka, = -2l£ ) (Q 2 )^(l ~x) + (Q 2 r (*) + (J - 5(1 - x) 



3^ , ^ m , . 3 1 - 2x + 
_ l v Q {x)+V 1 (x)-- + 



x 2 -x 3 



log(l — x) 



2\2 



(1 - X + x ) 

x(l — x) 



log(x)+0(e) 



(4.31) 
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Finally there are two gluon-gluon antennae 
1 1 



F 



9,99 



2(Q 



2\2 



8s 2 8t 2 8s 2 8u 2 8t 2 8u 2 12st I2su \2tu 
— + — + + + — + — + + — + 

t s 



4t 3 4u s 



4s 



u 



u 



u 



f — + — + — + 24s + 24t + 2Au + 0{e) 
su st tu 



1 /+2 
^9,qq (Q2)2 



^)+0( e ) 

S S ' 



Their integrated forms are given by: 

^,„ = -n2(e 2 )i(i-z) + (<? 2 )- 



11 2(1 - 2x + x 2 - x 3 ) , , , 2(1 -x + x 2 ) 2 , , 

+2P X x - — + ^ log 1 - x - log x 

ox x x(l — x) 



+0(e) 



2I« F (Q 2 )5(1 - x) + (Q 2 ) 



9,11 



2\-e 



b -5{l-x)+ l -V 1 {x) + ±- + 0{e) 



(4.32) 
(4.33) 

(x) 

(4.34) 
(4.35) 



5. Initial-initial configurations 

The last situation to be considered is when the two hard radiators are in the initial state. 
The subtraction terms necessary to account for singularities associated with these configu- 
rations are constructed in terms of initial-initial antennae. At NLO, one unresolved parton 
is emitted off these two radiators, as displayed in Figure 3. As before, more final state 
partons can be emitted at higher orders. 

The initial-initial configuration is slightly more involved than the previous two. Even 
though at NLO the integration of the antenna functions over the factorized phase space 
will turn out to be trivial, in order to guarantee this factorization, a very restricted kind 
of mappings will be allowed. In addition, to fulfill overall momentum conservation, both 
the hard radiators and all the spectator momenta, including non-colored particles, have to 
be remapped. This is done with a convenient generalization of the Lorentz transformation 
introduced in [3]. 

5.1 Subtraction terms for initial-initial configurations 

The NLO antenna subtraction term, to be convoluted with the appropriate parton distri- 
bution functions for the initial state partons, for a configuration with the two hard emitters 
in the initial state (partons i and k with momenta p\ and p 2 ) can be written as: 



dfiS,(ii) = ^ d$ m+ i(&i, . . . kj, k j+1 ,.. . , k m+1 ;p 1 ,p 2 ) 



m+l 



m+1 



^2 X ik,j(PliP2i k j) Mmih, • • • , kj-!, k j+1 , k m+1 ;X!P!,X 2 p 2 ) 



-Jm ^(^1> • • • i kj-li kj+!, • • • , &m+l, 



(5.1) 
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As mentioned, all the momenta in the arguments of the reduced matrix elements and the 
jet functions have been redefined. The two hard radiators are simply rescaled by factors 
x\ and X2 respectively. The spectator momenta are boosted by a Lorentz transformation 
onto the new set of momenta {ki, I / j}. As before, the mapping must satisfy overall 
momentum conservation and keep the mapped momenta in the mass shell. In this case, 
this turns out to severely restrict the possible mappings. 
We start from the (m + l)-parton phase space 

d^ m+1 (h,...,k m+ r,Pi,P2) = (2Tr) d 5 (^ Pl +p 2 -^k}j Y\[dki] (5.2) 



and insert 



and 



1 = J d d qd d q 5 (pi +p 2 - kj -q)S (xipi + X2P2 ~ q) , (5.3) 
1 = j n*(fc,-B(fc,,g,?))[dfc,], (5.4) 



where B is a Lorentz transformation that maps q onto q. We also insert 



1 = J dx\dx25(x\ — xi)S(x2 — X2) (5-5) 



with 

/ S\2 ~ Sj2 S 12 - Sij - S j2 



x 1 



S12 S12 - Slj 



1 



£2 = ( *u-°W-*V-°* )\ (5.6) 

V s 12 S12 ~ S j2 J 

These last two definitions guarantee the overall momentum conservation in the mapped 
momenta and the right soft and collinear behavior, they are derived in detail in Section 5.2 
below. Now we can integrate over the original momenta, k\ , I 7^ j by inverting the Lorentz 
transformation. The Jacobian factor associated with this integration is unity, as B is a 
proper Lorentz transformation. We also integrate over the auxiliary momenta q and q, to 
obtain 

d$ m+ i(fci, . . . , k m+1 ;pi,p 2 ) = d$ m (fci, . . . , kj-i, k j+1 , k m+ i; x±pi, X2P2) 

x8(xi — X\) 5{x2 — £2) [dkj] dxi dx2 ■ (5-7) 

At this point the phase space totally factorized into the convolution of an m particle phase 
space, involving only the redefined momenta, with the phase space of parton j. 

Inserting the factorized expression for the phase space measure in eq. (5.1), the sub- 
traction terms can be integrated over the antenna phase space. The integrated form of the 
subtraction terms must be, then, combined with the virtual and mass factorization terms 
to cancel the explicit poles in e. In the case of initial-initial subtraction terms, the antenna 
phase space is trivial: the two remaining Dirac delta functions can be combined with the 
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one particle phase space, such that there are no integrals left. We define the initial-initial 
integrated antenna functions as follows: 

Xikj(xi,x 2 ) = J [dk j ]x 1 x 2 5(xi -xi)6(x 2 -x 2 )X iktj (5.8) 

Substituting the one particle phase space, and carrying out the integrations over the Dirac 
delta functions, we have, 

Xikj(x 1 ,x 2 ) = (Q 2 y e T ^ _ ~j J(xi,x 2 )Q 2 X ikjj , (5.9) 

with Q 2 = q 2 = (pi +p 2 — kj) 2 . The Jacobian factor, J(x\,x 2 ) is given by 

„, \ x x x 2 (1 + x\ x 2 ) ( {1 + xi)(l + x 2 )\~ t . . 

J(X1 ' X2)= (x 1 + x 2 ) 2 {1 ~ Xl) (1 " X2) { (x 1 + X2 ) 2 ) ' (5 - 10) 

and the two-particle invariants are given by: 

x 1 (l-xl) x 2 (l-x\) 

Slj=S 12 ; , S j2 = -Si2 ; • (5.11) 

Xl + X 2 Xl + X 2 

The integrated subtraction term is, then, 

m+l j ^ m+1 J Xl X2 

Xd$ m (ki,.. .,kj-i,kj+i, . . .,k m+1 ;XiPi,X 2 P2) 

x \M m (h, k 31 ,k j+ i, k m+1 ;x 1 pl,x 2 p 2 )\ 2 

xJ^ih, kj-i, fcj+i, . . . , k m+1 ) , (5.12) 

where we have relabeled all k{ — > k{. The final step is to convolute this subtraction term 
with the parton distribution functions of the initial state particles. The integrated version 
of the subtraction pieces is then combined with the virtual and mass factorization terms 
to render a finite contribution when e — > 0. Recasting appropriately the convolutions, the 
integrated subtraction term is 



da s,(u) = 

m+l j 



±?^S m +iJ & J 6 Jfr x x 7 6 x 2 \x 2 J 



xC(e) X ikJ {xx,x 2 ) da B (txHx, &H 2 ) . (5.13) 
5.2 Phase space mapping 

By asking for momentum conservation and phase space factorization, we are severely con- 
straining the possible phase space mapping. The principal origin of this constraint is that 
the remapping of both initial state momenta can only be a rescaling, since any transversal 
component would spoil the phase space factorization. 

The two mapped initial state momenta must be of the form 

Pi = xxpx Pi = x 2 p 2 , (5.14) 
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so that 

q = Pl + P2 

is in the beam axis. Since the vector component of q = p\ + P2 — kj is in general not along 
the p\ — P2 axis we need to boost all the other momenta in order to restore momentum 
conservation. The transformation must map q onto q. As it must keep all the spectator 
momenta, which are arbitrary vectors, on mass-shell, it must belong to the Lorentz group. 
This transformation then fully determines the initial-initial phase space mapping, by fixing 
X12 in terms of the invariants. 

We consider a candidate Lorenz transformation A(q). It has to map the vector q into 
a vector A(q)q = q in the beam axis. From the result q of the transformation, one can read 
off xi and X2 using 

2(P1 +P2)q = Ol +X 2 )si2 
2(P1 - P2)q = (X2 - X±)s 12 



yielding 



Xl 



2p2q 

S12 
2 P i 



x 2 = (5.15) 

S12 

The two equations can be combined to give 

•512 - Slj - S 2 j 

x\x 2 = , 

S12 

which can also be derived from the on shell condition q 2 = q 2 . To ensure that the mapping 
has the right soft and collinear limits at NLO it is sufficient to impose A(q) = 1 for q in 
the beam axis. 

For the transformation A we take a boost Bj,(q) of appropriate parameter whose 
direction is transverse to the beam axis in the rest frame of Pi and P 2 . Objects defined in 
the rest frame of the new system are denoted by a *. This transformation clearly satisfies 
the requirement B^(q) = 1 for q in the beam axis, since then no boost is required to bring 
q in the beam axis. By construction, the longitudinal component of q in the rest frame of 
Pi and P2 is conserved, that is 

* (Pi-p 2 )g _ (Pi - m ,- 1fi . 

XlX 2 -\/S X\Xi\j s 

So that we have 

(x 2 - xi) = X2S2 > ~ X1S1 > (5.17) 
S12 

which gives the mapping 



Xi 



'-12 


- S 2 j 


js\2 


~ S lj ~ 


S 2 j 


S12 


~ S lj \ 




S12 




S12 




js\2 


~ S lj ~ 


S 2j 


«12 


- s 2j \ 




S12 





»» = i— \ J > ( 5 - 18 ) 
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which was used in (5.6) above. It yields the correct soft and collinear limits at NLO: 

1. j soft: xi — ► 1, X2 — ► 1. 

2. fej = zipi || pi: xi = (1 - zi), x 2 = 1. 

3. fej = z 2 p 2 || P2: = 1, x 2 = (1 - z 2 ). 

It should be pointed out the transformation is not unique. Possible transformations 
are however strongly constrained. If one requires a symmetrical treatment of x\ and x 2 , 
rotations are not allowed as transformation. To show that, we take pj to be transverse to 
the beam axis. Bringing q to the beam axis with a rotation will force us to choose to rotate 
q either towards the p\ or the p 2 side. This would favor either x\ or x 2 . The only way 
to bring q to the beam axis, without having to choose between x\ and x 2 is in this case a 
boost transverse to the beam axis. 

The extension of the phase space mapping to NNLO is trivial. In this case, four-parton 
antenna functions Xnjk require a mapping with two partons, j and k unresolved. The 
transformation B^(q) is unchanged, but now the vector q is given by q = p\ + p 2 — kj — kk- 
The momentum fractions in (5.6) are replaced by 



These two momentum fractions satisfy the following limits in double unresolved configu- 
rations: 

1. j and k soft: x\ — ► 1, x 2 — > 1, 

2. j soft and k k = z\p\ || p\: x\ — > 1 — z\, x 2 — > 1, 

3. fcj = zipi || pi and A: fc = z 2 p 2 \\ p 2 : x\ —> 1 - Z\, x 2 —> 1 - z 2 , 

4. % + k k = zipi || pi: xi -> 1 - zi, x 2 -»■ 1, 

and all the limits obtained from the ones above by exchange of p\ with p 2 and of kj with k k . 
The factorization of the phase space into an m-parton phase space and an antenna phase 
space goes along the same lines as for the NLO case. At NNLO, however, the integration 
of the antenna functions over this factorized phase space is no longer trivial. 

As in the initial-final case, we also require the NNLO mapping to turn into the NLO 
mapping (5.6) if only one parton becomes unresolved. In the limits where j becomes 
unresolved between i and k, we denote the parameters of the reduced NLO phase space 
mapping by x[ and x' 2 - We find: 




(5.19) 
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1. j becomes soft: 



i 

S12 ~ Sk2 S12 ~ S lk - S k2 \ 2 , 

XI -> [ = x 1: 

Sl2 Sl2 ~ Sl k J 

1 

r-2 - \ — — ! =x' 2 

Sl2 Si2 ~ S k2 



2. kj || fcfc, kj + h = K K : 



i 

2 



3. kj = zipi || p x : 



S\2 ~ SR2 S12 - SiK ~ S K 2 

xi -> =x 1 , 

S12 S12 - Sik 

Sl2 ~ SIR S12 ~ Sik - SR2 \ 5 , 

x 2 -> I = x 2 

«12 S12 - S K 2 J 



,'(1 - Zi)si2 - Sfc2 (1 - Zi)s 12 ~ (1 - Zl)Slfc - S fc2 \ - _ , . , 

X\ — > I — (1 — ZlJXj^ , 

«12 «12 - Slfe 



S12 - Slfc (1 - Zl)si2 - (1 - - S k2 



1 

2 



«12 (1 - Z 1 )S 12 ~ S k2 

All other single unresolved limits involving one radiator parton in the initial state follow by 
exchange of pi with p 2 or kj with k k . To subtract all unresolved limits of parton j between 
emitter partons i and k from Xu jk, one needs to subtract from it the product of an 
initial-final antenna function with an initial-initial antenna function Xijk • Xi^k- Analytic 
integration of these terms over both antenna phase spaces results in a double convolution 
in the rescaling variables for p\ and a single convolution in the rescaling variable for p 2 . 

At subleading colour, j can also become unresolved between i and I. In this case, we 
denote the reduced phase space mapping parameters by x'[ and x 2 . The limits read: 

1. j becomes soft: 

( S\2 ~ s k2 S12 - s lk - s k2 \ 2 „ 
xi -> =x 1} 

V s 12 Sl2 ~ Slk J 

( Sl2 ~ Si k Sl2 ~ Slk ~ Sk2 V // 

x 2 -> = x 2 

V s 12 S12 - S k 2 J 

2. fcj = zipi || pi: 

/(l - Zl)si2 - s k2 (1 - zi)su - (1 - zi)si fc - s fe2 \ 5 n w/ 
Xl -»• = (1 - Zi)Xi , 

V «12 «12 - Slk ) 

( Sl2 ~ Slk (1 - Zl)si2 - (1 - Zi)sife - S k2 \ 1 „ 
X2 r = S 2 ■ 

V -512 (1 - Zi)Si2 ~ S k2 J 
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3. kj = z 2 p 2 \\P2- 



S12- s k2 (1 - z 2 )s 12 - sik - (1 - z 2 )s k2 \ 2 _ „ 

' S12 (1 - 2 2 )S12 - Slfc / 1 ' 



f (1- z 2 )s 12 - s lk (I- z 2 )s 12 - si k - (1- z 2 )s k2 \ 2 n . „ 
x 2 -> = (1 - 2 2 )x 2 . 

V s 12 «12 - S k2 ) 

These single unresolved limits are subtracted from Xnj k by the product of two initial-initial 
antenna functions Xuj-Xjl^. Analytic integration of these terms over both antenna phase 
spaces results in two double convolutions in the rescaling variables for p\ and p 2 . 

5.3 NLO antenna functions 

The unintegrated antenna functions necessary to subtract all the singular configurations 
at NLO with two initial state hard radiators, can be obtained immediately from the corre- 
sponding initial-final ones quoted in section 4.3 by crossing. We have 

x9 kJ = S ikJ x9. k , (5.20) 

where 5i k j is an overall sign, which is —1 for D Q gg q and E®- ql and +1 for all other antennae. 

The Mandelstam variables of the unintegrated antennae in section 4.3 have to be 
replaced by s = (pi + p k ) 2 , t = (pi - pj) 2 , u = (p k - pj) 2 . 

Again, the splitting of antenna functions into different sub-antennae is different from 
the two configurations discussed above. For the initial-initial configurations there is no 
need to split the quark-gluon antenna D qg g as in all the singular limits it collapses to the 
same two-particle antenna. So D qg g is given by the crossing of (4.27). However, in this 
case, the antenna D gg q has to be split into two subantennae to separate the two collinear 
limits present in it. We find: 

D 99,q = D gi92,q + D 929l,q ' ( 5 - 21 ) 

such that D gig2 q contains only the singular configurations when the quark becomes collinear 
with gluon g 2 . Explicitly: 

n 1 ( s 2 4t 2 4u 2 3st 2t 3 3tu \ /nf . . 

As mentioned, the integration of the initial-initial antennae over the factorized phase space 
is trivial and only involves a proper treatment of the singularities when e — > 0. The 
integrated antennae read 



A = (O 2 



1 

2e Pqg 



-p q ° g ) {x 2 )5{l-x l ) 
1 1 — 2x 2 + 2xo / 1 — x 2 



H2 + ' ~T~M 2 tt^iw-^ ! 
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1 - 2x 2 + 2x\ xi(l + xix 2 ) {2xix\ + x 2 2 + (2x\ - 2x\ + l) x\) 



+ 



+ 



2(1 -xi) ' (xi+x 2 ) 3 (l-x 2 ) 
xix 2 (l + x\x 2 ) (2xi - x 2 (xf + 2x 2 xi — l)) " 



(xi + x 2 ) 3 



o 



r(l) 



2\-e 



(xf + x|) (x|xf + 2x 2 xf + xf + 3x 2 xi + 2xi + l) 



P u = 

19,9 



2(xi + l)(x 2 + l)(xi +x 2 ) 2 

((1 _ Xl )2 _ 2x 2 lQg ( IH-l) + (1 _ x 2) lQg (_JL^ S{1 _ X2) 
+ 2(1 -X X ) 

+^tt 2 5(1 - xi)<5(l - x 2 ) - 1(1 + x 2 )P (xi) + l£>o(xi)Po(x 2 ) 
+5(1 - xi)Pi(x 2 ) + (xi <-> x 2 ) 
-2(l«(Q 2 )) S(l-x 1 )S(l-x 2 ) 

,2\-e 



(Q 

2 



- ^rillMSQ. - x 2 ) - ^pg)(x 2 )5(l - xi) 



7T 



+— <5(1 - xi)<5(l - x 2 ) + 5(1 - x 2 )Pi(xi) + 5(1 - xi)Pi(x 2 ) 



+ l^xi + lo S V J 



1 — Xl 



1 + xi ( l-x x 
2 V xi + 1 



6(1 -x 2 ) 



f (-x 2 i + x 2 -2 + — )P (xi) 



1 + Xl 



^o(x 2 )+P (xi)P (x 2 ) 



(xix 2 - 1) (x|xf + (2x1 + ^2) x? + (2x| - xl + 4) xf ) (xix 2 + l) s 

xi (l - xf ) (xi + x 2 ) 3 (l - xf) 
(xix 2 - 1) (5xix 2 + 2x2) (xix 2 + l) 2 

xi (l - xf) (xi + x 2 ) 3 (l - xl) 
xl-x 2 + 2- ± 4 - xix 2 (x 2 + xi (xl + Xl + l) (x 2 + 1) - 3) 



1 — Xl 

2 



+ 



2xix 2 (l -xf) (l-xf) 



+ 



log(^) (4 - 4 + 2x2 - 1) log (2gf ) 



1 — x 2 

(Q 2 r 

2\-e 



X2 



5(1 -Xi) 



+(Q 
4 



lp(o) (x2 ) 5( i_ Xl ) 
l + l-2x 2 + 2x 2 



1 — x 2 



2 2 V 1 + x 2 

(2x 2 xf + x\ + xl (2x\ - 2x\ + !))(! + x 2 xi) 2 



5(1 -Xi) 



(1-x 2 ) (x 2 + xi) 4 



(5.23) 



(5.24) 



(5.25) 
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[kx\x\ + 3 (2xf + x i) x l + ( 4x i + 2x i + 1 ) x 2 + 3xfx 2 + xf) (x 2 xi + 1) 

2(x 2 + xi) 4 

(2xf + 2x^ - + 1) Pofo) 1 - 2x 2 + 2x| 1 



£ qq,q' ~ (Q 



£qq' ,q — (Q ) 



2(1 + xi) 2 2(1 -xO 

2 A - e ^l -2 (1 + xi x 2 ) 2 (1 - xf) 



2(xi +x 2 ) 4 
l p f)( Xl )5(l-x 2 ) 



+ (xi <-> x 2 ) , 



(5.26) 
(5.27) 



PoO^) (2 - 2xi + xf) x \ -2xi + 2 



2xi 2xi(l-x 2 ) 
x 2 (xix 2 + 1) ((xix 2 - l)xf + 2) 

xi(xi + x 2 ) 2 [x\ - l) 
(2xi + (2 - 2xi + xi) log (^1) - 2) 5(1 - x 2 ) 



2xi 

= -I^(Q 2 )5(1 - x 1 )5(l - x 2 ) + (Q 



(5.28) 



2\-e 



-lpg) (xi)( 5(l-x 2 ) 



+ ^vr 2 5(l -xi)<5(l -x 2 ) + f-xi + x 2 -2 + ^- ) 'A,:.n ) + •• .r 2 )Pj.(,r ;! ) 



+ 



1 + ( lo S " l) (*i ~4 + 2-2 - 1) (log (2^)) 



1 - x 2 



X2 



6(1 -xi) 



+ lp (xi)P (x 2 )- 2(i + J (i + x2) 



2(x 4 -xf + 2) (l + xix 2 ) 2 
+ (xi + x 2 ) 4 + 



x 2 - xi + 2 + -Q^fc - ± 2 (x 4 + x 2 + (xf + 3) xi + 2) (1 + xix 2 ) 5 



1 — x 2 



(1 -x%)(xi + x 2 y 



3 (x? + xi) [xl + x 2 ) (1 + xix 2 ) 2 { x lx\ + [x\ + l) 2 x\ + x 2 ) (1 + X1X 2 ) 



(xi+x 2 ) 4 
| 2xi (xf + 3) (l + xix 2 ) 2 | 



+ 



(1 + x x x 2 ) 



(x 2 + l)(xi + x 2 ) 



xi(xi + l)x 2 (x 2 + 1) 



(xi + x 2 ) 4 
+ [x\ <-> x 2 ) 



G°- 
y qq,9 



(Q ) 



2\-e 



7 gq,i 



= (Q 



2\-e 



(1 + xi x 2 )xf (1 - x 2 ) 2 (l + x 2 ) ; 

(Xi + X 2 ) 4 

- ^p£W(i-x 2 ) 



+ (xi <-> x 2 ) , 



(5.29) 
(5.30) 



2e 



(2(1 - xi) - (2 - 2xi + xj) log (2^))S(1-X 2 ) Vo(x 2 ) ( 2 - 2x x + xf) 



2xi 

2 - 2xi + xf + 2(1 + xix 2 ) (xf + 2x 2 xi + (xf - 2xf + 2) x 2 ) 



2xi 



2xi(l-x 2 ) 



2xi(xi + x 2 ) 3 (l - xl) 



(5.31) 
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6. Application of method at NLO 



To illustrate the application of the antenna subtraction method at NLO, we derive the 
antenna subtraction terms for two example reactions: deep inelastic (2+l)-jet produc- 
tion and hadronic vector-boson-plus-jet production. Several NLO calculations are already 
available in the literature both for the DIS jet production [3, 23] and for the vector boson 
production [24]. 

6.1 (2-|-l)-jet production in deep inelastic scattering 

The production of (2+1) jets in deep inelastic lepton-proton scattering can be described on 
the parton level by the scattering of a space-like virtual gauge boson and a parton, yielding 
a final state with two hard partons (with the extra jet coming from the proton remnant, 
not participating in the hard interaction). We limit ourselves to consider the scattering of 
transversely polarized virtual photons. The cross section can be written as 

da = j f E mo d ^ + fs(o ■ (e- 1 ) 

The partonic cross sections up to NLO are, in turn, given by, 

da q = d® 2 (k g ,k q ;p q ,q) |M°<J 2 jf\k g ,k q ) 

+d$ 2 (k g , k q ;p q , q) 2M(M^ gq M^ gq ) J 2 (2) (k g , k q ) 
+d^ 3 (k gi ,k g2 ,k q ; Pq ,q) \Mq >ggq \ J 2 ^ (k gi , k g2 , k q ) 
+d<&3(k qi ,k q2 ,k q ;p q ,q) \M qqqq \ J 2 \k q , k q , k q ) 

+ ^2d<f> 3 (k q ,k q/ ,k q r, Pg ,q) \M° qqlql \ J^\k q , k q >,k q >) , (6.2) 
q'+q 

da g = d<S> 2 (k q ,k q ;p g ,q) \M^ qq \ 2 jf\k q ,k q ) 

+d$ 2 (k q , k q ; p q , q) 2»(M„ 1 ) t M° ) J 2 (2) (k q , k q ) 



gm am 

+d<S> 3 (k g , k q , k q ;p q , q) \Ml gqq \ 2 J { 2 \k g , kg, k q ) . (6.3) 



For the sake of brevity, the momentum arguments on the matrix elements were omit- 
ted. The first index on the matrix elements indicates the incoming parton, the remaining 
ones the outgoing partons. Subtractions for infrared real radiation singularities must be 
performed only on the three-parton final states. The matrix elements for the real contri- 
butions can be expressed in terms of colour-ordered three-parton and four-parton antenna 
functions. They read as follows: 

WlJ = e 2 q N 2<q A° 3 (l q ,3 g ,2 q ) (6.4) 

I 2 _ 2 1 



\ M q,ggq\ = <^ q 



NA$(l q , 3g,4 g , 2g) + N A u 4 (l q , 4 g ,3 g ,2 q ) 

(6.5) 



A\(l q ,3 g ,<± g ,2q) 
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M 



o 

<?><?<?? I 



+ \ M q,qq'q'\ 



e 2 q N 3 , q 



N F B 4 (l q ,3 q ,,4q>A) 



--jy(C®(l q , 3 g , 4g, 2q) + C2(l q , 3 q , 2g, Aq) 

+ C° A (4q, 2q, l q , 3 q ) + C°(2 q , 4 q , l q , 3 q )) 

+ ^e 2 q , N 3 , q BZ(3q,,lq,2q,4q,) , 
q' 



(6.6) 
(6.7) 



M, 



o 



g,gqq\ 



3,9 



N A°(l q , 3 g , 4 g , 2q) + TV Al(l q , 4 g , 3 g , 2 5 ) 



iV 



^4(l<?;3 9 ,4 9 ,2 g ) 



where 



and 



/\r 2 - 1 

iV n , g = 4™ ( 5 2 )- 2 ^^ 2(1 - e)Q 2 , 



iV n , 9 = 4™ ( 5 2 r" 2 2Q 2 . 



(6.8) 



(6.9) 



(6.10) 



All the antenna functions used here can be obtained from explicit expressions for the final- 
final antennae X 3 and X® in [14] by crossing in each case the particle denoted by hi. 

The subtraction terms for both quark and gluon initiated processes are a combination 
of final-final and initial- final subtractions. We split the quark-induced contributions into 
three terms: quark-plus-two-gluon final states at leading and subleading colour, da^ A and 
dff^ and quark-quark-antiquark final state da^ B . Identical-only quark contributions to 
the matrix elements, involving the antenna C®, have no single collinear limits so they do 
not need to be subtracted. Gluon-induced contributions can also be split into leading and 
subleading colour, der*^ and da^. The explicit expressions for subtraction terms are 



given by 



d °q,A = el N 3,qd<f>3(k gi ,k g2 ,k q ;p q ,q) 



N 



x { D l 9l92 A lGq4 2 \K G ,k g ) + Dl g2q Al GQ 4 2) (K G ,K Q )) , (6.11) 



= -e 2 q N 3: qd<f> 3 (k gi ,k g2 , k q ;p q ,q) 



2N 



{ A i 91 q A Q, 92 Q 4 2 \k 9 „K Q ) +Al g2q A° Q ^ Q 4 2 \k gi ,K Q )) , (6.12) 



&°q,B = N 3,qd^3{k q ,k q/ ,k q r,p q , 



+ 4 { E lq>q> A Q,Gq A^^G, k q ) + A^q jf\K G , k q )) 



Y,4> E lq> A G,Q>q>4 2) (K Q ,,k q ,) 



(6.13) 
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^2 e l N 3,gd$3(k g , k q , kq-;p g , q) N 



g,A 



D g,gg ^Q,Gg 'h ( K G, k q ) + D ggq -A GqQ - J 2 (k q , Kq) 

4 N 3t gd$ 3 (kg, kg, kg^g, ?) ^ 



(6.14) 



x {^ AQ g,gg A Q,gQ J 2 + ^O.sQ ^9> ^q) 



(6.15) 



In the above, Xj t j^ and X^jj denote three-parton antenna functions with momenta /, J 
obtained from a phase space mapping. The combination of these subtraction terms with 
the real matrix elements containing three partons in the final state is finite in all soft and 
collinear limits and can be integrated numerically over the three-particle phase space. 

On the other hand, the analytical integration of the subtraction terms over the factor- 
ized phase space can be carried out using the results of Section 4.3. For the poles of the 
integrated terms, we obtain: 



Ka + ^ q , A + Kb 



= -2 



a s 
2vr 



+N F fl£> >F (t) + daf(p,g) 



(i) 



da% A + da s , 

9< A g,A 





f dx 1 


2vr J 


x e 




f dxl 


2vr J 


x e 


a s 




2 — 


N (l 


2vr 





(6.16) 



+Nf [i qg ]At) + ^l F (u))]d*Z( P , q ) 



qg,F^ 





f dx 1 


2vr J 


x e 


a s 


r dxi 


2^ J 


x e 



(6.17) 

where Cf = (N 2 — 1)/(2N), Tp = \ and the Born cross sections are given by 

daf(p,q) = d$ 2 (k g ,k q ;p,q) \M° qm \ 2 J?\k g ,k q ), (6.18) 

(6.19) 



daf (p,q) = d<S> 2 (k q ,k q ;p,q) \M^ J^'(k q ,k q ). 



q,gq\ 

|2 T (2), 



The poles contained in the operators 1^ match exactly the ones appearing, with opposite 
sign, in the interference of the renormalized one loop amplitudes with the Born ones. The 
remaining poles correspond to the mass factorization contributions. Thus, combining the 
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integrated subtraction terms with the virtual contributions and the mass factorization 
counterterms, we obtain a finite contribution, free on any poles in e that can be integrated 
over the two partons phase space. 

6.2 Vector-boson-plus-jet production in hadronic colliders 

The second example we will consider is the production of a vector boson V (V = 7, Z, 
W^) plus a hadronic jet in a hadronic collision. This process is mediated by the scattering 
of two partons into the vector boson and one hard parton. The cross section is given by 



da= r |i |2 r ^ [hMi)f _ (6) + f _ Mi)f (&)] d& 
j a 6 ^ 

+ E [4(6)4(6) + 4(6)4(6)] da, i9j . 

+ E K4(6) + 4(6)) / 3 (6) + / 9 (6) (4(6) + 4(6))] da, 

i 

+/ 3 (6)/ 9 (6)da 99 }. 



(6.20) 



Again we express the partonic cross sections in terms of color ordered antennae. We first 
write: 

2 



d 4<?j = d< ^2{kg,q;p q ,Pq) 



.\I 



_ 



J^) 



+a$ 2 (k 9 , q ;p q ,p q )2M(M 1 ^ jtg M^ j j4 1 \k 9 ) 



+d<Z> 3 (k gi ,kg 2 ,q;p q ,p q ) 
+ ^2d<Z>2(k q ,kq,q;p q ,pq) 



M°- 



j[ \k gi , k 92 ) 



k,i 



M, 



o _ 



Ml 



o 



k,l 



d 4s = Z~2,d<S> 2 (k q ,q;p q ,p g ) 



Ml 



o 



qiQj,qkqi 



4 1] (k q ) 



J\ {kqz i k qi ) , 
J\ ^ ( k q3 ' kq 4 ) , 



(6.21) 
(6.22) 



+ £ d<& 2 (fc„ ? ; p„p 9 ) 25R(M^ i9 .M ? ° iffi? .) jf > (fc, 



d<r, 



35 



+ Xl d$ 3(^ k g ,q;p q ,p g ) 
j 

^2d<f> 3 (k q ,k q ,q;p gi ,Pg 2 ) 



M° 

q t g,qjg 



Ml 



o 

gg,q%qj 



j[ \k q , kg) 

j[ ^ (k q , k q ) , 



(6.23) 



where we have omitted again the momentum arguments of the matrix elements. The matrix 
element for the partonic process ab — > cdV is given by M^^ and the momentum of the 
vector boson, appearing in the phase space measure is denoted by q. 
The real contributions are given by 

= |Cy 2 {vf + a^N^AlQ^lg^) (6.24) 



Ml 



o 



q t qj,g 
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\Ci 3 \ 2 K 2 + a 2 )N lm A° 3 (j q ,i g ,i q ) 



\Cij\ 2 (^ + a?)^o, w 



N A\ (j q , l g , 2 g , i 5 ) + N A° 4 (j q ,2 g ,l g , i q ) 



Aiijq, Igj 2 S , i g ) 



(6.25) 



(6.26) 



N 2 ,qq {5 k l \Cij\ 2 (vf + af) feg, ig, *g) 

+<% |Cjw| 2 (vl + al)Bl(k q ,j q ,iq,lq) 
+Sji \C ik \ 2 (vj + af) B2(k q ,j q ,lq,iq) 
+$ik \Cji\ 2 {v 2 + a 2 ) B%(j q ,k q ,i q ,lq) 



Vi Vk &ly{jq, kq, l q , i q ) + a% Ofc B^ A {jq, k q , l q , i, 
+5i k 5jl 2 ^(CaCjj) (vi Vj B^ v (kq, j q , l q , iq) + Oj Gtj B2 A (k q ,j q , lq, Iqfj 

2 

+ ^ <5»fc<5 fcZ \Cij\ 2 (vj + af ) C£(j g , k q , l q , i q ) 
2 

+ N S i3 d 3l \ C ik\ 2 {V 2 + a 2 )Cl{kq,] q ,lq,iq) 

2 

+^ fytffci \Cij\ 2 (vj + af) C^(i 5 , lq, k q ,j q ) 



2 i 

+j-f5ijdik \C jk \ 2 (v 2 + a 2 )C%(l q ,i q ,k q ,j q )j 



(6.27) 



7V 2 , 



?9 



2 



{^7 ICjfel 2 {vf + a 2 ) Bl(k q ,l q ,j q ,i q ) 
+Sik \Cji\ 2 (v 2 + a 2 )Bl(l q ,k q ,i q ,j q ) 

+5 ik 5ji 2 ^(CjjC^) (v, Bly{k q , lg,jg, iq) + a, aj B^ A (k q , l q ,j q , iqj^J I 
Ski 



+N 2 ,qq [Sij (1 - Ski) + 

x{s jk \C U \ 2 (v 2 + a 2 )B2(lq,kq,j q ,i q ) 
\C jk \ 2 (v 2 + a 2 )Bl(k q ,l q ,i q ,j q ) 

+6j k 6u2$l(CiiCj j ) (vi Vj Bly{lq, k q ,jq, iq) + Oj dj B^ A (l q , kg,jq, iq)^J 

2 

-Jf 5 3l 5 3k \ C ij\ 2 ( v i + ° 2 ) C2(«5, fc <?> U 

2 

—jySiiSji \Cik\ 2 (vj + a 2 )Cl{k q ,l q ,i q ,j q ) 
2 

-jjSikSjk \Cu\ 2 {v 2 + a 2 )Cl(l q ,k q ,i q ,j q ) 

~ 6 ik Sii |Cy 2 («? + af) Cj&, fc„ I,)} , (6.28) 
|Cyf K 2 + af) ^ [iV^(i ? ,i ff ,2 ff ,y + AT^(j„2 s ,i ff J ? -) 



M, 



o 



9l92,qiqj 



\C i3 \ 2 (v 2 + a 2 )N 2 , gg \N A° 4 (i q ,i g ,2 g ,j q ) + N Al(i q ,2 g ,i g ,jg) 



— A^iq, lg, 2g,j q/ 



where 



Nn,qq — ^n,qq 



4ira v (g ) 



N n , qg = 4na v {g 2 ) n ' 2 -Q'- 



2\n-2 



N 2 



N 2 



(l-e)Q 2 



N n , qg = 4 Tray 
with the coupling constants given by 



2\n-2 



N 2 — 1 1 — e 



a- 



a 



a w 



Q, 



az 



G F M 2 z ^/2 
64tt 



47T ' " 47T 

The vector and axial couplings of the quarks to the vector bosons are 

7 2 7 1 

vl = V v d =--, 



(6.29) 

(6.30) 

(6.31) 
(6.32) 
(6.33) 

(6.34) 



0. 



l-^sm 2 9 w , v$ = -1 + ^sm 2 6 w , 



(6.35) 



V U — V d — ! a il — a d — ~ ■ 

The flavor mixing matrices, Cjj are given by <5jj in the case of 7 and Z production and by 
the CKM matrix in case of W production. Finally, the colour ordered antenna functions 
appearing in eqs. (6.24) to (6.30) can be obtained from explicit expressions for the final- 
final antennae X® and X® in [14] by crossing the particles denoted with hats. 

The subtraction terms for this process involve initial-initial and initial-final antennae 
as there are at most two partons in the final state, final-final antennae are not needed. Only 
antennae A® and B® contain singular configurations and, thus, require to be subtracted. 
We find the following subtraction terms, classified according to the partonic reaction they 
must be combined with, 

1 



da 



QiQ j, 99 



\Cij\ (Vi +a i )-N 2 ,qqd^ 3 (kg 1 ,kg a ,q;p q ,pq) 



aa qiqj,q k qi 



-Jr [^4,fl,,o a J i W ^ (6-36) 

^2 N 2 , qq d^ 3 (k q ,k q ,q;pq,pq) 
{S kl \C tJ f (v 2 + of) \ (El q ,-, A%^ G J? {K G ) 



k,l 

X 



- 28 - 



da. 



s 

QiQj,1k<ll 



+5 3 i \C lk \ 2 (vf + ahE -^ A% iG a k J?\K Qk ) 
+5 lk \C^ (v* + a*)E^ qk A% jG ^\K Ql )}, 

QiPquPqj) 



(6.37) 



k,l 



{ ( x - t ) ^ IQfc|2 + <o, Qk 4 1] (K Qk) 

+ (l - I) ** l^| 2 (y? + a?)*^ < GiQi J?W) 

+ (<% (i - ««) + ^ i^i 2 ^ 2 + ° 2 ) A lo, Ql 4 1] (k Qi ) 

(6.38) 

X] \ C H I 2 K 2 + ° 2 ) d$ 3 (A; ? , kg,q;p q ,p g ) 
j 



qi9,qj9 



+ d U 9 a Iq j , g 4 1) (k g ) 



1 r 

"N 



9i92,qiqj 



1^' I 2 + ° 2 ) ^2-99 d$ 3(fe g , %, q;p 9 i,Pg 2 ) 



(6.39) 



iV 2 -l 



Antennae of the form Xij t x and Xjj^k correspond to antennae where the momenta of the 
particles denoted with capital letters are obtained by initial-final and initial-initial phase 
space mappings respectively. 

The integrated form of the subtraction terms can be obtained inmediately using the 
results in Sections 4.3 and 5.3. The singular pieces of these terms are then given by 



da s 



qiqj,gg 



= -2- 



d<7, 



qiqj,qkqi 



aa qiqj,qkqi 
^°qi9,q]9 



= -2 



'2vr 

a s f dx 1 
'2tt 



a s 
2— 
2vr 

a s 

2^ 

a s 

2vr 

a s 



N (l qg \t) + I$(u)) -— do5 ft (p 9 ,p ff ) 



^-C F p^{x) (da* q .(xp q ,p q )+daf iq: .(p q ,xp q )) + 0(e°) ,(6.41) 



r(!) 



/ T^ F * (a:) (<^(P 9 .^)+d^p(PJ.^))+0(6°), (6-42) 
/ ^C^JJO*) (d^(^,xp 9 ) + d^(p 9 -,x^))+0( e °), (6.43) 
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a s f dx 1 m B 

a s f dx 1 
2it J x e 



a, f dx 1 



) (x)+Np^ F (x)) da* g (p q ,x Pg ) 

2 / _ T^gW^d^.fe.x^ + O^), (6.44) 

j 

d ^,M, = -£ / ^7 2TF * (X) £ ( d ^^) + d^(xp g -,p 9 )) + O(e ) , 

i 

(6.45) 

where the Born cross sections are given by 



d< ^9 = 5^ d *2(fc„9;P„P S 



_ 

9i<?j>9 



^1 ^(^9) ! 



M° 

qig,qj 



(6.46) 
(6.47) 



Again, the poles in the operators 1^ are canceled by the virtual contributions whereas 
the ones associated to the Altarelli-Parisi kernels drop out once combined with the mass 
factorization counterterms. 



7. Conclusions and Outlook 

In this paper, we have generalized the antenna subtraction method for the calculation of 
higher order QCD corrections to exclusive collider observables to situations with partons 
in the initial state. 

The basic ingredients to the subtraction terms, the antenna functions, can be ob- 
tained from the known final-state antenna functions by simple crossing. We derived the 
factorization of an multi-parton phase space into an antenna phase space (required for the 
analytic integration of the subtraction terms) and a reduced phase space of lower multiplic- 
ity, for antennae with one or two hard radiator partons in the initial state (initial-final and 
initial- initial antennae). Explicit phase space factorization and parameterization formulae 
were presented for NLO and NNLO calculations. We derived all integrated initial-final 
and initial-initial antennae relevant at NLO, and demonstrated their application on two 
example calculations. 

A major advantage of the antenna subtraction method is its straightforward extension 
to NNLO calculations. Our results are a significant step towards NNLO calculations of 
hadron collider observables. Using the phase space factorizations presented here, NNLO 
subtraction terms for jet production observables at hadron colliders can be constructed 
from known building blocks. Their analytic integration over the antenna phase spaces 
relevant to NNLO calculations is still an outstanding issue. It is however anticipated that 
usage of techniques similar to those applied for the integration of the final-final antennae 
will help to perform these integrals in a systematic and efficient way. 

First applications of the method presented here, once the corresponding NNLO in- 
tegrated antennae are available, could be NNLO calculations of two-jet production or 
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vector-boson-plus-jet production at hadron colliders, and of two-plus-one-jet production 
in deep inelastic scattering. Further extensions of the method could include radiation off 
massive particles, thus allowing the NNLO calculation of top quark pair production at 
hadron colliders. 

Another important extension of subtraction methods is the combination with parton 
shower algorithms [25], thus allowing for a full partonic event generation to NLO accuracy. 
This task was fully accomplished so far only for one QCD subtraction method [4]. While 
fixed order NLO calculations are independent of the subtraction method used, there can 
be a residual dependence on the method in matched NLO-plus-parton-shower calculations, 
since unintegrated and integrated subtraction terms are treated differently in the parton 
shower. With the formulation of the antenna subtraction method for initial state radiation 
presented here, it will become possible to construct antenna-based parton showers for 
hadronic collisions. 
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